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We introduce a three-state model for a single DNA chain under tension that distinguishes between 
B-DNA, S-DNA, and M (molten or denatured) segments and at the same time correctly accounts 
for the entropy of molten loops, characterized by the exponent c in the asymptotic expression S ~ 
— clnn for the entropy of a loop of length n. Force extension curves are derived exactly employing 
a generalized Poland-Scheraga approach and compared to experimental data. Simultaneous fitting 
to force-extension data at room temperature and to the denaturation phase transition at zero force 
is possible and allows to establish a global phase diagram in the force-temperature plane. Under 
a stretching force, the effects of the stacking energy, entering as a domain-wall energy between 
paired and unpaired bases, and the loop entropy are separated. Therefore we can estimate the loop 
exponent c independently from the precise value of the stacking energy. The fitted value for c is 
small, suggesting that nicks dominate the experimental force extension traces of natural DNA. 



I. INTRODUCTION 



DNA continuously stays in focus of polymer scientists due to its unique mechanical and structural properties. In 
particular the possibility to trigger phase transformations in this one-dimensional system has intrigued theorists from 
different areas [l|. In fact, the thermal denaturation or melting transition of DNA was shown to correspond to a 
true phase transition, brought about by a logarithmic contribution to the configurational entropy of molten loops or 
bubbles, S ~ —clnn, as a function of the loop size n 0|. The value of the exponent c is crucial since it determines the 
resulting transition characteristics. For c = 3/2, the value for a phantom chain without self-avoidance, the transition 
is continuous, while self-avoidance increases c slightly beyond the threshold c = 2 above which the transition becomes 
discontinuous (HQ. A distinct mechanism for transforming DNA involves the application of an extensional force. For 
forces around F w 65 pN DNA displays a highly cooperative transition and its contour length increases by a factor of 
roughly 1.7 to 2.1 over a narrow force range 4 6]. These experiments sparked a still ongoing debate on whether this 
over-stretching transition produces a distinct DNA state, named S-DNA, or merely the denatured state under external 
tension. According to the first view S-DNA is a highly stretched state with paired bases but disrupted base stacking 
In the other view the over-stretched state consists of two non-interacting strands (l4l - [r7| . Evidence for the 
existence of a distinc t S-state co mes from theoretical models @, [1] , molecular dynamics simulations fl0l - fl3T | and from 
AFM experiments of iRief et al.l fl8l - l2"oT ] where in addition to the over-stretching transition a second weak transition at 
forces between 150 pN and 300 pN is discerned, which has been interpreted as a force induced melting of the S-state. 
The critical force of both transitions depends on the actual sequence [l9[ and the salt concentration [2l|, but the 
interpretation of the second transition is complicated by the occurrence of pronounced hysteresis effects that depend 
on various parameters such as pulling velocity, salt concentration, or presence of co-solutes such as cisplatin jl9l [20|. 
On the other hand, support for the view according to which S-DNA is not a distinct state comes from theoretical 
models [3 US], simulation s I22I a nd recent experiments by Shokri et al. [IH and van Mameren et al. 17 1. 



Apart from simulations |1CH13| |. existing theoretical works that grapple with experimental force traces or DNA 
melting fall into three categories with increasing computational complexity, for reviews see refs. [23U25I ] . In the first 
group arc Ising-likc models for DNA under tension which give excellent fitting of the over-stretching transition but 
by construction cannot yield the denaturing transition [l5l. |26| - |28| . The work of Marko [2!| is similar but employs a 
continuous axial strain variable. In the second group are models that include a log arithmic entropy contribution of 
molten loops in the spirit of the classical model by Poland and Scheraga @ @, IH, |30h3^ | . This gives rise to effectively 
long-ranged interactions between base pairs and thus to a true phase transition. The third group consists of models 
which explicitly consider two strands [35j-l39| . Those models thereby account for the configurational entropy of loops 
- at the cost of considerable calculational efforts - and correspond to loop exponents c = d/2 in the absence of 
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self-avoidance effects, where the dimensionality of the model is d = 3 for rcf. [35U37j and 1 for rcf. [H, |39|. All these 
above mentioned works consider only two different base states (paired versus impaired) and thus do not allow to 
distinguish between B-DNA, S-DNA and denatured bases. Recently, three-state models were introduced that yield 
very good fits of experimental force traces at ambient temperatures. However, in previous analytic treatments of such 
three-state models 0, l40| , the loop entropy was neglected and therefore the temperature- induced denaturation in the 
absence of force cannot be properly obtained, while the loop entropy was included in a simulation study where most 
attention was given to dynamic effects @ . 

In this paper we combine the Poland-Schcraga formalism with a three-state transfer matrix approach which enables 
us to include three distinct local base pairing states and at the same time to correctly account for long-ranged 
interactions due to the configurational entropy of molten DNA bubbles. Our approach thus allows for a consistent 
description of thermal denaturation and the force induced BS-transition within one framework and yields the global 
phase diagram in the force-temperature plane. We derive a closed form expression for the partition function of 
three-state DNA under tension. This allows to systematically investigate the full parameter range characterizing the 
three states and the DNA response to temperature and external force. In our model we allow for the existence of 
S-DNA but stress that the actual occurrence of S-DNA is governed by the model parameters. By assuming such a 
general point of view, our work is able to shed new light on the question of the existence of S-DNA. The extensible 
worm-like chain model is employed for the stretching response of each state. The loop exponent is found to have quite 
drastic effects on the force extension curve. For realistic parameters for the stacking energy, the experimental force 
extension curves are fitted best for small loop exponents < c < 1, hinting that the DNA in the experiments contained 
nicks. Loop exponents c > 1, which give rise to a genuine phase transition, are not compatible with experimental 
force-distance curves. Under external force, the effects of stacking energy and loop exponent arc largely decoupled, 
since the stacking energy only determines the coopcrativity of the BS-transition while the loop exponent influences 
the second SM-transition found at higher forces. This allows to disentangle these two parameters, in contrast to the 
denaturation transition at zero force where the effects of these two parameters are essentially convoluted. The precise 
value of c is important also from a practical point of view, as it impacts the kinetics of DNA melting [jj |4l| , which is 
omnipresent in biological and bio-technological processes. 



II. THREE-STATE MODEL 



Double-stranded DNA is modeled as a one-dimensional chain with bases or segments that can be in three different 
states, namely paired and in the native B-state, in the paired stretched S-state, or in the molten M-state. The free 
energy of a region of n segments in the same state reads 

Ei{n,F) =n- gi {F)-5 lM -k B T\nn- c , (1) 

with i = B, S, M. The force F dependent contribution 

9i{F) = 9° + 9 r ctch (F) + 9 y hC (F) (2) 

is split into three parts, g® is a constant that accounts for the base pairing as well as the difference of reference 
states of the worm-like stretching energy, cf. supporting material eq. (??). The stacking energy of neighboring bases 
in the same state is absorbed into g®, too, so that the stacking energy will appear explicitly only as an interfacial 
energy Vij between two regions which are in a different state. The second term g^ trctch = —F 2 li/(2-Ki) takes into 
account stretching along the contour with li and Ki the segmental contour length and the elastic stretch modulus. 
Finally, gf /hC (F) is the free energy of a worm- like chain (WLC) in the Gibbs ensemble (constant force F), based on the 
heuristic relation between force F and projected extension x F i WLC (a;) • ^/keT = (l — x/(n^)) /A+x/(nli) — I /4 
where £j is the persistence length and n the number of segments. The Gibbs free energy n ■ gJ^ LC (F) of a stretch of n 
segments is extensive in n and follows via integration, see supporting material section ??. We note that this is only 
valid if the persistence length is smaller than the contour length of a region, £j < nli, which is a plausible assumption 
because of the high domain wall energies. Likewise, the decoupling of the free energy into contour stretching elasticity 
and worm-like chain elasticity is only approximate p3l |44| but quite accurate for our parameter values [HI : For small 
force WLC bending fluctuations dominate and the contour extensibility is negligible, while contour stretching sets in 
only when the WLC is almost completely straightened out. The last term in eq. |T]) is the logarithmic configurational 
entropy of a molten loop (i = M), characterized by an exponent c [3l I30L l46j . see supporting material section ??. The 
exponent is c = 3/2 for an ideal polymer [I?} and 2.1 for a self avoiding loop with two attached helices [U, |48[. If 
the DNA loop contains a nick the exponent is reduced to c = for an ideal polymer and 0.092 for a self avoiding 
polymer [3[. We consider the simple case c = 0, where transfer matrix methods can be used to yield results in the 
canonical ensemble with a fixed number of segments N Q , as well as the case of finite c where we introduce a modified 
Poland-Schcraga method to obtain results in the grand canonical ensemble. 
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III. PARTITION FUNCTION 
A. Modified Poland-Scheraga approach for c / 

The molecule is viewed as an alternating sequence of different regions each characterized by grand canonical partition 
functions. Various techniques for going back to the canonical ensemble are discussed below. The canonical partition 
function of a stretch of n segments all in state i = B, S, or M is 

Qi(n) = exp(-0Ei(ri)) , (3) 

where j3 = (keT) -1 is the inverse thermal energy. The grand canonical partition functions are defined as Z{ = 
S^Li ^ n Qi( n ) with A = exp(/3/i) the fugacity and fx the chemical potential. The grand canonical partition function 
of the whole DNA chain which contains an arbitrary number of consecutive B, S, and M stretches reads 

oo 

Z = v T • (Mp S Vp S ) fc M PS • v = v T • (1 - Mp S Vp S )- 1 M PS • v, (4) 

fc=0 

with the matrices in this Poland-Scheraga approach given by 

<Z B \ / e-^BM^ 

Mp S = | Z s , V PS = c-? VsB e~P VsM \ , v = ( 1 ) (5) 
Z M J y c -?v MB e -/3V MS o 

and where 1 is the unity matrix. The energies V^- are the interfacial energies to have neighboring segments in different 
states and are dominated by unfavorable base pair un-stacking. The diagonal elements of Vps are zero which ensures 
that two neighboring regions are not of the same type and thus prevents double counting. The explicit form of Z 
is given in the supporting material, see eq. (??). The partition function in eq. (U) is general and useful for testing 
arbitrary models for the three DNA states as given by the different Zi. This approach is also easily generalized to 
higher numbers of different states. Using the parameterization eq. ([!} for vanishing loop exponent c = the partition 
functions of the different regions arc given by 

°° Ac^^ 1 

Zi = £ \ n Qi{n) = 1 _ Xe _ 0gi , for Ae-"* < 1, (6) 

n—1 

i = B, S, M. Insertion into eq. (j4]) yields 
a 1 \ + a 2 \ 2 +a 3 A 3 

which is a rational function of the fugacity A, whose coefficients - determined by eqs. ((4]) and (|6]) - are smooth 
functions of the force F and the temperature T. For c ^ the partition function of a molten stretch is modified to 




oo oo 

£A"Q M (n) = £A»(e^)"l 

n—1 n—1 

Li c (Ae~^ M ) , for Ae~' 39M < f (8) 



where Li c (z) = Y^=i z " n c f° r ^ < f is the polylogarithm [4!| and exhibits a branch point at z = 1. The functional 
form of the grand canonical partition function for c ^ reads 

= & A + M 2 + 6 2 Li c (A/A b ) + MLi c (A/A b ) + o 4 A 2 Li e (A/A b ) 
c# ° 6 5 + 6 6 A + 6 7 A 2 +o 8 ALi c (A/Ab) + & 9 A 2 Li c (A/Ab) ' U 

where Ab = c^ 9M denotes the position of the branch point and the coefficients 6j, determined by eqs. © and ©, 
are smooth functions of F and T. 

The grand canonical ensemble where N, the total number of segments fluctuates, does not properly describe a DNA 
chain of fixed length. We therefore have to investigate the back-transformation into the canonical ensemble where the 
number of segments iV is fixed. For the back-transformation there are three options: 
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1. Calculus of residues route: 

The grand-canonical partition function Z(X) = Y2n=i ^ N Q{N) can be viewed as a Laurent series, the coefficients 
of which are the canonical partition functions Q(N) determined exactly by 

ow-ssif^- (10) 

The contour C = Xoe 2mt , < t < 1, is a circle in the complex plane around the origin with all singularities of Z(X) 
lying outside. This complex contour integral can be evaluated using calculus of residues (Hoj which becomes technically 
involved for large N and thus limits the practical relevance of this route. 



2. Legendre transformation route: 

The canonical Gibbs free energy 

g(N) = -k B T In Q(N) (11) 
and the grand potential 

$(/i) = -k B TlnZ(A), (12) 
are related via a Legendre transformation 

g{N) = $(n(N)) + N'ft(N). (13) 
The chemical potential /iasa function of the segment number N is obtained by inverting the relation 



Let us briefly review the origin of eqs. (fl3| and fj 14[) in the present context. Changing the integration variable in 
cq. (|10p from A to \i = ln(A)//3, the complex path integral can be transformed into 

Q(N) = [ e-P^-^^&n w e -PH»*pW)-pNu sp (N) ^ ^ 



Jc 

with the contour C = fXo + 2irit/f3, < t < 1 and /io = keTlnAo- The integral in eq. tjl 5|) has been approximated by 
the method of steepest descent, where the contour C is deformed such that it passes through the saddle point /z sp |50| | 
determined by equation (fT4")) . If $ features singularities deformation of the contour C requires extra care. In the 
present case the presence of a pole A p = exp(/3/Lt p ) of Z(X) produces no problem as /x sp < jjL v holds, meaning that the 
deformed contour does not enclose the pole singularity. This is different for the branch point singularity /ib where we 
will encounter the case /ib < A*sp f° r large c > 2. 



3. Dominating singularity route: 

For large systems, i.e. N S> 1, one approximately has — In Q(N) ~ iVlnAd, where the dominant singularity 
Ad = exp(/3/xa) is the singularity (in the general case a pole or a branch point) of Z(X) which has the smallest 
modulus. One thus finds 

G{N) =k B TNlnX d . (16) 

This easily follows from eq. (fTU)) : In the limit of JV" S> 1 the integral can be approximated by expanding Z(X) around 
Ad and deforming C to a Hankel contour which encircles Ad [El[. F° r the case where N(fi) oc (/id — /i) _Q , a > 0, this can 
be understood also in the context of a Legendre transform. As TV = —d$>/djj, one has $ oc (/id — /i) _Q+1 and therefore 
the first term of cq. (H3J) scales like $(/x(AT)) cx N 1 ' 1 / 01 . Thus, the second term Nfi(N) cx N^ d - TV 1 " 1 /" « 7\r^ d i s 
dominant. Since the saddle point behaves as /i sp = fi(N) — > /id for N — > oo, it follows that the dominating singularity 
expression cq. (|16p equals the Legendre transform cq. (|13[) in the thermodynamic limit N ^ oo. 
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B. Transfer matrix approach for c = 

For c = only interactions between nearest neighbors are present and straight transfer matrix techniques are 
applicable. We introduce a spin variable i n for each segment which can have the values i„ = B, S, M. The energetics 
are given by the Hamiltonian 

N N-l 

H(ix,i 2 ,...,iN) = y^,9i„ + ( 17 ) 

n—l n—1 

where gi n and Vi n i n+1 are the previously introduced parameters for the segment and interfacial free energies. The 
canonical partition function of the molecule can be written as 

Q(N)= e-^ il ' i2 '-' iw >=v T -T ,v - 1 M T M-v, (18) 

n »« 

where we introduced the transfer matrix T = MtmVtm and 
/ e -f>9B \ / 1 

M TM = c"^ 3 , V TM = ( e^ VsB e -W SM | v = ( 1 | . (19) 

\ c -PawJ 

Q(N) is calculated readily by diagonalizing T 

Q(N) = v T • UD^^U^Mtm • v = v t . D^- 1 ■ v r = £ v^xf ~ l . (20) 

where D = U _1 TU is a diagonal matrix with eigenvalues Xi, the columns of U are the right eigenvectors of T, 
•vf = v T • U and v r = U~ Mtm • v. By virtue of the Perron- Frobenius theorem one eigenvalue x max is larger than 
the other eigenvalues and thus the free energy is in the thermodynamic limit dominated by x max and reads 

Q = -k B T\nQ(N)^-k B TN\nx max . (21) 

The transfer- matrix eigenvalues Xi and the poles A p of Z c= q eq. ([7]) are related via Xi = 1/A p . As expected, the free 
energies from the transfer matrix approach eq. (|21j) and from the Poland-Scheraga approach for c = and using the 
dominating singularity approximation eqs. (|16j) are identical in the limit N — > oo. Clearly, for c ^ the modified 
Poland-Scheraga approach yields new physics that deviates from the transfer-matrix results. 

Although not pursued in this paper, the transfer matrix approach allows to calculate correlators. For example the 
probability pu{k, m) of a denatured region with k consecutive molten base pairs starting at base m is given by 

p M (k, to) = Q(N)' 1 ■ v T ■ T m - 2 (P B T + P s T)(P M T) fc (P B T + PsT)T iV_m_fc_1 MTM ■ v . (22) 

The Pi-matrices, which project a segment onto a certain state, are defined as 

/0 0\ 

and P M = . (23) 
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IV. FORCE EXTENSION CURVES 

The central quantity is Q (F, T, N) , the Gibbs free energy of a DNA chain with N base pairs, subject to a force F and 
temperature T. From Q (F, T, N) , obtained via the Legendre transform, eq. (fT3)) , the dominating singularity, eq. p^|) , 
or the exact transfer matrix partition function, eq. (|20|) . we can calculate observablcs by performing appropriate 
derivatives. The number of segments in state i — B, S, M is obtained by 



(24) 

T,N,F 



The force extension curve is readily calculated via 



E ||= E N^x^^+F-k/^ (25) 



T.N „ % dF 



=B.S,M J i=B,S,M 



where xJ^ hC (F) is the stretching response of a worm-like chain and given explicitly in the supplementary information 
cq. (??). 
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Figure 1: Comparison of force extension curves obtained by different methods for c = 0. The curve obtained via the exact 
transfer matrix calculation eq. (|20|) is already for N = 2 accurately reproduced by the approximate Legendre transformation 
eq. (|13p . The dominating singularity method eqs. (|16[1 or - equivalently - (|2ip is strictly valid in the thermodynamic limit 
but agrees with the Legendre transform already for a modest value of N = 10. The units of the abscissa is extension per base 
pair (bp). Parameters for A-DNA in the absence of DDP are used, see supporting material section ??. 

A. Vanishing loop entropy, c = 

In this section we compare the prediction for vanishing loop exponent c = to experimental data and obtain 
estimates of the various parameters. We also demonstrate the equivalence of the grand canonical and canonical 
ensembles even for small chain length N . 

In order to reduce the number of free fitting parameters we extract as many reasonable values from literature as 
possible. For the helical rise, the stretch modulus, and persistence length of B-DNA we use Zb = 3.4 A, kb = 1 nN, 
and £b = 48 nm [52[. For the M-state, which is essentially single stranded DNA (ssDNA), ab initio calculations yield 
lu = 7.1 A and km = 2 • 9.4 nN (53|, where km is valid for small forces F < 400 pN and the factor 2 accounts for the 
presence of two ssDNA strands. Our value for the stretch modulus is considerably larger than previous experimental 
fit estimates [3, HH which might be related to the fact that experimental estimates depend crucially on the model 
used to account for conformational fluctuation effects; however, the actual value of km is of minor importance for the 
stretching response, see supporting material. The persistence length of ssDNA is given by £m ~ 3 nm (55j . It turns 
out that the quality of the fit as well as the values of the other fit parameters are not very sensitive to the exact value 
of the persistence length £s and the stretch modulus ks of the S-state as long as 10 nm < £s ^ 50 nm and «s is of the 
order of km, see supporting material section ??. Therefore we set £s = 25 nm, which is an intermediate value between 
the persistence lengths of ssDNA and B-DNA, and «s = k m = 2 • 9.4 nN Q. The segment length of the S-state Is 
will be a fit parameter. 

The chemical potentials gf, i = B, S, M, account for the free energy of base pairing and, since we set the interaction 
energies between neighboring segments of the same type to zero, Vn = 0, also for the free energy gain due to base pair 
stacking [56|. They also correct for the fact that the reference state of the three different WLCs, which is x = in 
the Hclmholtz ensemble (constant extension x), cf. eq. (??), is not the same as contour and persistence lengths differ 
for B-, S-, and molten M-DNA. We choose g B = and treat g§ and g M as fitting parameters. 

Each of these parameters controls a distinct feature of the force-extension curve: The chemical potentials g® 
determine the critical forces, the segment lengths U affect the maximal extensions of each state and the off-diagonal 
Vij control the cooperativity of the transitions, see supporting material section ?? for an illustration and section ?? 
for a summary of all parameter values. 

1. Force extension curve 

In fig. [T] force extension curves based on three different levels of approximation are compared, using the same 
parameters that we extracted from DNA stretching data as will be detailed below. It turns out that the force 
extension curve obtained via the Legendre transformation route eq. (|13|) (dashed line) is a very good approximation 
of the results obtained from the exact transfer matrix results eq. (|2"0"|) (dotted line) already for N = 2. For N = 10 
and larger virtually no differences between these two approaches are detectable. The deviations from the dominating 
singularity route eq. (fT6]l (solid line), which gives a result independent of N, are somewhat larger. But one sees that 



Transfer matrix, N — 2 

■ Legendre transform, N = 2 

■ Legendre transform, iV = 10 

■ Dominating singularity 
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Figure 2: Bottom panel: Force extension curve of double-stranded A-DNA with and without DDP. Symbols denote experimental 
data [2C|. lines are fits with the three-state model for c = 0. The main difference between the two curves is the lack of 
cooperativity in the BS-transition in the presence of DDP which we take into account by choosing vanishing interaction 
energies Vy = 0, i,j = B, S,M. Top panel: Fraction Ni/N of segments in the different states as follows from eq. (|24[) in the 
absence of DDP. 



already the Legendre transform for N = 10 (dash-dotted line) matches the dominating singularity result very closely. 
Therefore the use of the dominating singularity, eq. (|16p , or the largest transfer-matrix eigenvalue, eq. (|2ip , is a very 
good approximation already for oligo-nucleotides and will be used in the rest of this work. 

In fig. [5] experimental stretching curves of A-DNA with and without DDP (cisplatin) are presented [2(J. When 
B-DNA is converted into S-DNA or M-DNA the base stacking is interrupted which give s rise to an interfacial energy 
between B and S as well as between B and M of the order of the stacking energy (Iollll|. For untreated DNA, we use 
the value Vbs = Vsb = Vsm = Vmb = 1-2 • 10~ 20 J and show in the supplement that variations down to 0.8 • 10~ 20 J 
do not change the resulting curves much. Vsm is presumably small as the stabilizing stacking interactions are already 
disrupted [3] , we thus set Vsm = Vms = for the fits in fig. [2] - but we will come back to the issue of a non-zero Vsm 
later on. Cisplatin is thought to disrupt the stacking interaction between successive base pairs and thereby to reduce 
the cooperativity of the BS-transition. This fact we incorporate by setting all interfacial energies to zero, Vy = 0, 
for DDP treated DNA. The three remaining undetermined parameters (Is, <?g, g M ) have distinct consequences on the 
force-extension curve. The segment length Is and the chemical potential g® determine the position of the BS-plateau 
with respect to the polymer extension and applied force, respectively, while the chemical potential g^ controls the 
force at which the second transition appears. Fitting to experimental data is thus straightforward and yields for 
untreated A-DNA l s = 6.1 A, g° = 1.6 • 10~ 20 J, g M = 2.4 • 10~ 20 J and for A-DNA in the presence of DDP (cisplatin) 
k = 6.0 A, g% = 1.2 • 10~ 20 J, g M = 2.8 • 10~ 20 J, see fig. H We also fit the number of monomers N and allow for 
an overall shift along the x-axis. The main difference between the two stretching curves is the cooperativity of the 
BS-transition, which is controlled by the interfacial energies Vbs and Vbm- Note that, although the over-stretching 
transition is quite sharp for DNA without DDP, it is not a phase transition in the strict statistical mechanics sense. 
A true phase transition arises only for c > 1, as will be shown in the next section. The fitted value of g M is about 
two times larger than typical binding energies [56| for pure DNA. As a possible explanation, we note that the force 
extension curve of DNA without DDP exhibits pronounced hysteresis (especially at higher force) which will increase 
the apparent binding energy due to dissipation effects [4(| ■ Any statements as to the stability of S-DNA based on our 
fitting procedures are thus tentative. However, such complications are apparently absent in the presence of DDP Q 
which rules out kinetic effects as the reason for our relatively high fit values of <7 M and the stability of S-DNA. Cisplatin 
most likely stabilizes base pairs due to cross-linking and thus shifts the subtle balance between B-, S-, and M-DNA. 
Therefore, the relative stability of B-, S-, and M-DNA is sensitively influenced by co-solute effects. We note that even 
with c = a good fit of the data is possible. In the top panel of fig. [2] we show the fraction of segments in B-, S-, 
and M-states for untreated A-DNA. There is a balanced distribution of bases in all three states across the full force 
range, in agreement with previous results Q. 



B. Non- vanishing loop entropy, c/0 



We now turn to non-zero loop exponents c ^ and in specific try to estimate c from the experimental stretching 
data. The partition function Z c ^q in eq. ([9]) exhibits two types of singularities. First, simple poles at A = A p , which 
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are the zeros of the denominator of cq. (j9)) and which are determined as the roots of the equation 



65 + b 6 X + M 2 
b s \ + b 9 \ 2 



Li c (A/A b ) . 



(26) 



Second, a branch point that occurs at 



A = A b = c 



Pgu 



(27) 



The singularity with the smallest modulus is the dominant one |l|, [5l|, and we define the critical force F c as the 
force where both equations, cq. (|26|) and (|27|) . hold simultaneously. For c < 1 the dominant singularity is always 
given by the pole A p and thus no phase transition is possible. For 1 < c < 2 a continuous phase transition occurs. 
By expanding eq. (p6|) around F c one can show that all derivatives of the free energy up to order n are continuous, 
where n <G N is defined as the largest integer with n < (c — 1) _1 [3]. For instance, c = 3/2 leads to a kink in the 
force extension curve. For c = 1.2 this leads to a kink in x'"(F). If c > 2 the transition becomes first order and 
the force extension curve exhibits a discontinuity at F = F c . In fig. we plot force extension curves for different 
values of the loop exponent c with all other parameter fixed at the values fitted for untreated DNA. It is seen that 
finite c leads to changes of the force extension curves only at rather elevated forces. In order to see whether a finite c 
improves the comparison with the experiment and whether it is possible to extract the value of c from the data, we 
in fig. [3}d compare the untreated DNA data with a few different model calculations for which we keep the parameters 
^Si ffs> ^M' ^BSi Vbm fixed at the values used for the fit with c = in fig. [5J Allowing for finite c but fixing a zero 
domain wall energy between S- and M-regions, Vsm = 0, leads to an optimal exponent c = 0.6 and slightly improves 
the fit to the data which show the onset of a plateau at a force of about 100 pN. The same effect, however, can be 
produced by fixing c = and allowing for a finite Vsm, which yields the optimal value of Vsm = 1.1 ■ 10~ 21 J. Finally, 
fixing Vsm = 1-1 • 10~ 21 J and optimizing c yields in this case c = 0.3 and perfect agreement with the experimental 
data. However, the significance of this improvement is not high, as the experimental data are quite noisy and possibly 
plagued by kinetic effects. What the various curves illustrate quite clearly, however, is that a non-zero exponent c 
leads to modifications of the stretching curves that are similar to the effects of a non-vanishing domain wall energy 
Vsm- Although Vsm should be considerably smaller than Vbm or Vbs, a finite value of Vsm = 1-1 ■ 10~ 21 J as found 
in the fit is reasonable and cannot be ruled out on general grounds. The maximal value of c is obtained for vanishing 
Vsm and amounts to about c = 0.6. A value of c = 2.1, which would be expected based on the entropy of internal 
DNA loops 0, on the other hand does not seem compatible with the experimental data, as follows from fig. [3^,. This 
might have to do with the presence of nicks. Nicks in the DNA drastically change the topology of loops and result in 
a reduced loop exponent which is c = for an ideal polymer and 0.092 for a self avoiding polymer Therefore, the 
low value of c we extract from experimental data might be a signature of nicked DNA. Additional effects such as salt 
or co-solute binding to loops are also important. Therefore c can be viewed as a heuristic parameter accounting for 
such non-universal effects as well. We note in passing that c only slightly affects the BS-transition, as seen in fig. [3^,. 
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Figure 3: Various force-extension curves of the three-state model with fit parameters for A-DNA without DDP. a) Lower 
panel: Force extension curves for different values of the loop exponent c, showing no phase transition (c < 1), a continuous 
(1 < c < 2), or a first order phase transition (c > 2). The critical forces are denoted by open circles. The inset is a magnification 
of the region around the transition. Upper panel: Fraction of bases in the three states for c = 3/2. The critical transition, 
above which all bases are in the molten M-state, is discerned as a kink in the curves, b) Comparison of experimental data 
(circles) and theory for c / 0. The curve for c = and Vsm = 0, already shown in fig. [H is obtained by fitting Is, gs, gm to 
the experimental data, the values of which are kept fixed for all curves shown. The curve for Vsm = and c = 0.6 results by 
fitting c and slightly improves the fit quality. The curve c = and Vsm = 1-1 ■ 10 -21 J is obtained by fitting Vsm- The curve 
for Vsm = 1-1 • 10 -21 J and c = 0.3 is obtained by fitting c and keeping Vsm fixed. The inset shows the first derivative of x(F) 
and illustrates that increasing c leads to a growing asymmetry around the transition region, c) Temperature dependence of the 
force extension curves. Increasing temperature leads to a decrease of the BS-plateau force. In the presence of a true denaturing 
transition, i.e. for c > 1, the critical force F c decreases with increasing temperature and for F > F c the force extension curve 
follows a pure WLC behavior. 



V. FINITE TEMPERATURE EFFECTS 



The temperature dependence of all parameters is chosen such that the force extension curves at T = 20 °C that 
were discussed up to now remain unchanged. The persistence lengths arc modeled as £,i(T) = ■ (T/293 K)~ The 
S-state free energy is split into enthalpic and entropic parts as q §(T) = rs(hs — Tss) where we use hs = 7.14 ■ 10~ 20 J 
and ss = 1.88- 10~ 22 J/K from Clausen-Schaumann et al. [lj|. The correction factor rs = 0.98 accounts for slight 
differences in the experimental setups and is determined such that g$(T = 20 °C) equals the previously fitted value. 
The molten state energy g M (T) = (/im — Tsm) is also chosen such that g M (T = 20 °C) agrees with the previous fit 
value and that the resulting denaturing temperature in the absence of force agrees with experimental data. Assuming 
a melting temperature of T c = 348 K for A-DNA H3, we obtain g M (T) = 1.5 • 10~ 19 J - T-4.2 ■ 10" 22 J/K for c = 
and <7m(T) = 1-6 ' 10" 19 J — T ■ 4.6 ■ 10 -22 J/K for c = 3/2. In fig. we plot a few representative stretching curves 
for different temperatures. It is seen that increasing temperature lowers the BS-plateau and makes this transition 
less cooperative. Differences between c = and c = 3/2 are only observed at elevated forces, where for c = 3/2 one 
encounters a singularity characterized by a kink in the extension curves. 

For c > 1 the critical force F c is defined as the force where the pole and the branch point coincide and eqs. (|2l)|) 
and (|27[) are simultaneously satisfied. The phase boundary in the force-temperature plane is thus defined by 

_ 6 5 + M b + Mg =c(c); 

8 Ab + OgA b 

where hi and Ab depend on T and F and £(c) = Li c (1) is the Riemann zeta function. The phase boundary F C (T) 
for c = 3/2 is shown in fig. [4] and agrees qualitatively with the experimental data of ref. [3]. For exponents c < 1 
no true phase transition exists. We therefore define crossover forces as the force at which half of the segments are in 
the molten state or in the B-state, i.e. Nm/N = 1/2 or Nb/N = 1/2. In fig. |4]we show these lines for both c = 
and c = 3/2. Note that the parameters for the c = case have been adjusted so that Nm/N = 1/2 at T = 348 K 
and F = 0. The broken lines on which Nb/N = 1/2 for c = and c = 3/2 are virtually the same, showing again 
that loop entropy is irrelevant for the BS-transition. The S-state is populated in the area between the Nb/N =1/2 
and Nm/N = 1/2 lines, which for T > 330 K almost coalesce, meaning that at elevated temperatures the S-state is 
largely irrelevant. Force- induced re-entrance at constant temperature is found in agreement with previous two-state 
models [30l [35| . Re-entrance at constant force as found for a Gaussian model [30l | is not reproduced, in agreement 
with results for a non-extensible chain [35| . 

As we have shown so far, a non-zero loop exponent c only slightly improves the fit of the experimental stretching 
data and the optimal value found is less than unity. This at first sight seems at conflict with recent theoretical 
work that argued that a loop exponent larger than c = 2 is needed in order to produce denaturation curves (at zero 
force) that resemble experimental curves in terms of the steepness or cooperativity of the transition 0] . To look into 
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Figure 4: The solid line shows the critical force F c for c = 3/2, at which a singularity occurs according to eq. (|28[) . Phase 
boundaries for c = (thick lines) and c = 3/2 (thin lines) are defined by Nm/N = 0.5 (dotted) and Nb/N = 0.5 (dashed). 
The melting temperature T c is denoted by the dot. The insets show the behavior of the phase boundaries near the melting 
temperature, F oc y/T c — T. Parameters for A-DNA without DDP are used. 



Figure 5: Relative fraction of segments in the B-state, Nb/N, as a function of temperature for different loop exponents 
c = 0, 1.5, 2.1 and for finite BS interfacial energy Vbs = Vbm = 1.2 ■ 10 -20 J (bold lines) and for Vbs = Vbm = (thin lines). 
Circles indicate the positions of the phase transition. For all curves parameters for A-DNA without DDP have been used, 
g M (T) = 1.5- 10~ 19 J-T-4.2- 10~ 22 J/K for c = and £&(T) = 1.6- HT 19 J - T ■ 4.6 • 10" 22 J/K for c> 0. 

this issue, we plot in fig. [5] the fraction of native base pairs, Nb/N as a function of temperature for zero force and 
different parameters. As soon as the domain-wall energies Vbs and Vbm are finite, the transition is quite abrupt, 
even for vanishing exponent c. Therefore even loop exponents c < 1 yield melting curves which are consistent with 
experimental data, where melting occurs over a range of the order of 10 K [57l. [H^]. 



The fact that the domain- wall energy due to the disruption of base pair stacking, Vbs and Vbm, and the loop 
entropy embodied in the exponent c, give rise to similar trends and sharpen the melting transition has been realized 
and discussed before 0, HU . The present three-state model and the simultaneous description of experimental data 
where the denaturation is induced by application of force and by temperature allows to disentangle the influence of 
these two important effects. By the application of a force, the de-stacking and the loop formation occur subsequently, 
allowing to fit both parameters separately. As our main finding, we see that for a finite domain- wall energy Vbs = Vbm, 
the additional influence of the loop exponent on the force stretching curves and the denaturation curves is small. In 
fact, the optimal value for c turns out to be of the order of c ss 0.3 — 0.6, even if we choose a vanishing value Vsm = 0. 
This estimate for c is smaller than previous estimates. One reason for this might be nicks in the DNA. So it would be 
highly desirable to redo stretching experiments with un-nicked DNA from which the value of c under tension could be 
determined. A second transition at high forces of about F ~ 200 pN which is seen in the experimental data used in 
the paper, inevitably leads via the fitting within our three-state model to an intermediate S-DNA state. But we stress 




Temperature T / K 



VI. CONCLUSIONS 
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that the occurrence of such an intermediate S-statc depends on the fine-tuning of all model parameters involved, 
which suggests that in experiments the S-state stability should also sensitively depend on the precise conditions. 
One drawback of the current model is that sequence effects arc not taken into account. This means that our fitted 
parameters have to be interpreted as coarse-grained quantities which average over sequence disorder. Calculations 
with explicit sequences have been done for short DNA strands but should in the future be doable for longer DNA as 
well. 
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The supporting material shows the calculation of the of Gibbs free energy gT rLC (F) of one worm- 
like chain segment subject to a force F. Additionally the explicit form of the grand canonical 
partition function Z is given. The loop exponent of an ideal polymer is derived and effects of self- 
, avoidance on the loop exponent are discussed. We illustrate that different model parameters are 

■ linked to certain features of the force extension curve. Tables with all parameters used in the main 
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GIBBS FREE ENERGY OF A WORM-LIKE CHAIN 



> 

in . 

To calculate the worm- like chain Gibbs free energy g™ hC (F), we first calculate the Hclmholtz free energy 
■ of a worm-like chain (WLC) in the fixed extension ensemble (Hclmholtz ensemble). After that we Legendre transform 
(x) in order to obtain (F), which is the thermodynamic potential of a WLC in the fixed force ensemble 

£S1 ■ (Gibbs ensemble). As we will see, Gf /LC (F) is extensive with respect to the chain's contour length nli, li is the length 
of one segment, and thus we can define the Gibbs free energy per segment gi VLC {F) = Gf /LC (F)/n. We introduce 
the dimcnsionless variables 



£i = 4> Fi = F-fo, #wlc = h wlc.^|1 ) GWLC^WLC.^ z = B,S,M, (SI) 
nli nli nli 

where is the persistence length of DNA in state i = B, S, M. The interpolation formula of Marko and Siggia [l[ for 



the force extension curve of the WLC 

= 4(1^ + (S2) 
is used to calculate the free energy in the Helmholtz ensemble 

S? LC (5) = £ F(x')dx' + HZ LC = + 4T C ■ (S3) 
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H^ LC is a free energy offset that accounts for the fact that even in the absence of an external force the free energies 
of chains consisting of B-, S-, or M-segments are not the same. In fact, this constant is not easy to calculate but can 
be dropped in the following since it can be adsorbed into the free energy contributions cf. eq. (??) in the main 
text. 



Now, we invert eq. (|S2 



4Fj _ i 
3 



1/3 



V 




gi {F) ~-n-w' G {Fm ~ m — Wi wm - 1) — m * ■ Xi {m * 



with ii(Fi) from eq. (|S4 



(S4) 



and perform a Legendre transformation to obtain the Gibbs free energy of a WLC in the Gibbs ensemble 

Gf LC (£) = H^ixtim - Ft ■ . (S5) 

Eqs. JSTJ and (J35) lead to 



(S6) 



B. EXPLICIT FORM OF THE GRAND CANONICAL PARTITION FUNCTION 



In the main text the expression 

oo 

Z = v T ' (Mp S Vp S ) fe M PS • v = v T • (1 - Mp S Vp S )- 1 M PS • v 



fc=0 



has been derived for the grand canonical partition function. The matrices are given by 



'Z B 
Mpg = | Z s 1, V PS = I c 
Z M Ve" 



q-PVbs e -0V BM - 

n p-pvs*, 





/3Vmb g-/3V M s o 



|1 



(S7) 



(S8) 



1 is the unity matrix. Evaluating eq. (|S7|I yields the grand canonical partition function of the three-state model 



Z B +Z S + Z M + {e~ pVBB +c^* aB ) Z B Z S + (e-^" SM + c^" MS ) Z S Z M + (c 
+ ( c ~ ,3(VBS+ysM) + e -P( v su+v M B) _|_ e -p(v M B+v BS ) + e -P{v B M+v us ) + e -p(v MS +v S B) _)_ e -/3(y S B+v BM ) 
e -P{v BS +VsB) _ e -^(VsM+v M s) _ e ~' 9(v ' MB+v ' BM)N )ZB^s-2 : M 



Z B Z S + (e"^ SM + c 



pv MS \ ZsZm + ( e -pv M B + e -Pv BM \ ZmZb 



I _ e -P(VBs+v S B) ZB2s _ c -P(v sm +v ms ) ZsZm e -p(v M B+VB M )z M z B 

-\ -1 

- ^ c ' i3 '- Vbs+Vsm+Vmb ^ + e~ l3 ^ VBM+VMS+VsB) ^jZ B ZsZi i 
which depends on the grand canonical partition functions Zi of the different regions. 



(S9) 



C. ORIGIN OF THE LOGARITHMIC LOOP ENTROPY 



In the simplest case a polymer can be described as an ideal random walk with step length b in d = 3 dimensions. 
When speaking of polymers b is also called Kuhn length. Let us first consider a one-dimensional random walk. The 
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probability to be after N steps at point x = nb is given by the binomial distribution 

p( T ) - n (N+n)/2 . ( , _ UN-n)/2 _ ( N \. n (N+n)/2 . /•, _ \(N-n)/2 (a 1() \ 

F[X} - ((N + n)/2)\((N-n)/2)\ P U P) ~{(N + n)/2j P (1 P) 

where p = 1/2 is the probability to move to the right and (1 — p) the probability to move to the left. By virtue of the 
central limit theorem the binomial distribution can be approximated by the normal distribution if N is large 

F(i)-^^e"w . (Sll) 

Now, we consider a random walk in three dimensions, where there are N/3 steps in each spatial direction. Since the 
steps in each of the three directions are independent from each other the probability distribution for being at point 
R = (x, y, z) after N steps is 

P(R) = P(x)P(y)P(z) = p^ e-l^ . (S12) 

The probability of a closed random walk, i. e. a loop, is the probability to return to the origin after N steps 

P ^ - (nmIw • (S13) 

The entropy difference between an ideal polymer forming a loop and an unconstrained polymer is therefore given by 
AS loop {N) = k B lnP(0) - k B lniV- c , (S14) 

where we introduced the loop exponent c = c\deaX = d/2 = 3/2 

If one considers now self-avoiding polymers the loop exponent increases, csaw = dv £± 1.76 with v ~ 0.588 in d = 3 
dimensions Q. Hence the entropic penalty for closing a loop increases. However, helices emerging from the loop 
increase c further. In the asymptotic limit of long helical sections, renormalization group predicts q = dv + 07 — las 
for a loop with I emerging helices 0, 0] where m = el(2 - l)/16 + e 2 l(l - 2) (81 - 21)/512 + C(e 3 ) in an e = 4 - d 
expansion. One obtains c\ = 2.06 for terminal loops present in hairpin structures and c-i = 2.14 for internal loops. 

The same analysis holds if one considers different topologies, namely peeling off one strand from the other starting 
from a nick. In this case one can show that c = for an ideal polymer and c = 0.092 for a self avoiding polymer 0,0]. 

Application of force to loops also alters the value of the loop exponent, which has been shown by Hanke et al. 0. 



D. EFFECT OF PARAMETERS ON FEATURES OF THE FORCE-EXTENSION CURVE 

Different parameters of the model are linked to certain features of the force extension curve. The chemical potentials 
g® determine the critical forces (fig. 1ST]) , the segment lengths U affect the maximal extensions (fig. |S2|) of each state 
and the off-diagonal Vij determine the cooperativity of the transition (fig. IS3[) . 
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Figure SI: The chemical potentials gf influence the transition forces. The solid line shows the force extension curve for the 
parameters used in the main text for DNA without DDP. The dashed and the dotted line show stretching curves with different 
chemical potentials of the S-state <?g while all other parameters are kept constant. Lowering g2, leads to a decreasing plateau 
force of the BS-transition as the free energy difference between B- and S-states (gs(F) — <?b(-F)) decreases, see eq. (??) in the 
main text. At the same time the force of the SM-transition increases as the free energy difference between S- and M-states 
(9m(F) - gs(F)) increases. 




Figure S2: The segment lengths U determine the extensions between which the transitions occur. The solid line shows the 
force extension curve for the parameters used in the main text for DNA without DDP. The other curves have been obtained 
by by varying the segment length of the S-state Is while keeping all other parameters constant. The effect of changing Is is 
twofold. First, increasing Is elongates the BS-plateau. As a side effect a larger Is stabilizes the S-state as its gs rLC (F) is 
lowered. This can also be concluded from eq. (|S6[) . Thus, increasing Is shifts the equilibrium towards the S-state, which is the 
same as decreasing gs , cf. fig. IS1I 
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Figure S3: The interfacial energies Vij determine the cooperativity of the transitions. The solid line shows the force extension 
curve for the parameters used in the main text for DNA without DDP. The other force extension curves have been obtained 
by varying Vbs or Vsm while keeping all other parameters constant. Lowering Vbs (dotted, dashed, and dash-dotted lines) 
smoothes the BS-transition and renders it less cooperative, which has been observed for DNA with DDP, cf. fig. ?? in the main 
text and the experiments by Krautbauer et al. [f|. Vsm has similar effects on the melting transition: Increasing Vsm from its 
original value to 1.2 • 10~ 20 J (gray curve) sharpens the denaturation transition. 
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Figure S4: The influence of the stretch modulus on force extension curves is quite moderate. The solid line shows the force 
extension curve for the parameters used in the main text for DNA without DDP. The dashed and the dotted line show the 
data with km = 40 nN and km = 9.4 nN, respectively, while keeping all other parameters unchanged. One observes that 
changing km by a factor 4 results in deviations of the extension by just a few percent even at extreme forces F = 400 pN. The 
dash-dotted line shows a force extension trace where ks = km = 1.6 nN, and is = Im = 0.59 nm have been used, which are the 
typical experimental values for two parallel ssDNA 0]. With this choice of parameters the force cannot discriminate between 
the S- and M-state anymore - the two states are degenerate. Therefore we have to set the chemical potentials equal, fit them 
again and obtain = $ M = 1.6- 10 -20 J. Consequently the shoulder indicating a denaturation transition does not appear in 
the force extension curve in contrast to the other curves shown. 
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E. TABLES WITH MODEL PARAMETERS 



Parameter 


bymbol 


A-DNA A-DNA + 


JDP Reference 


segment length 


7 

'B 


0.34 


0.34 


Wenner et al. |8| 


(nmj 


Is 


0.61 


0.60 


fitted 




hi 


0.71 


0.71 


Hugel et al. [9] 


persistence length 


SB 


48 


48 


Wenner et al. \o\ 


(nmj 


>- 


25 


25 


assumed 






6 


■_> 

6 


Murpny et al. |1U] 


stretch modulus 




1 n 
l.U 


l.U 


Wenner et al. ^ 


l niN J 


Kg 


O O A 


o o a 
z • y.4 


nugei ci ai. y 




km 


2-9.4 


2-9.4 


TTno-ol r»t ftl rQl 
nugei ei di. ^7 


chemical potential 


9°b 


0.0 


0.0 


defined 


(io- 20 J) 


9°S 


1.6 


1.2 


fitted 




9°M 


2.4 


2.8 


fitted 


interfacial energy 


Vbbj T^sSi Vmm 


0.0 


0.0 


incorporated in g,° 


(io- 20 J) 


^BS 


1.2* 


0.0 


assumed, disrupt stacked base pairs [11] 




Vbm 


1.2* 


0.0 


assumed, disrupt stacked base pairs [11] 




Vsm 


0.0 


0.0 





Table SI: Parameters for stretching curves at temperature T — 293 K. 



Parameter 


Symbol 


A-DNA 


Reference 




persistence length 




48/(T/293 K) 






(nm) 


&(r) 


25/(T/293 K) 
3/(T/293 K) 






chemical potential 




= hi - Tsi 






enthalpy 


ft,B 


0.0 






(io- 20 J) 


hs 


7.0 


Clausen-Schaumann et al. 


[12] . slightly rescaled 




h M 


15 


fit for c = 






hyi 


16 


fit for c = 3/2 




entropy 


•SB 


0.0 






(10- 22 J/K) 


•ss 


1.8 


Clausen-Schaumann et al. 


[121. slightly rescaled 




SM 


4.2 


fit for c = 






sm 


4.6 


fit for c = 3/2 





Table S2: Temperature dependence of parameters for DNA without DDP. The enthalpy hyi and entropy sm of the denatured 
state are determined by fixing the denaturation temperature T c for zero force at T c = 348 K. 
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